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a  =  Dimensionless  semi-minor  axis  of  ellipsoid  (cl./D")  , 
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'a  =  Semi-minor  axis  of  ellipsoid,  speed  of  sound 
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f  =  Dimensionless  function 
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Relative  Mach  number  of  drop  C^a."*^//^; 


M        =   Shock  wave  Mach  number 
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m  =  Mass  of  drop 
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T  =  Dimensionless  time  for  shock  wave  to  pass 
over  drop 

T  =  Dimensionless  time  for  boundary  layer  separation 
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2  =  Condition  downstream  of  shock  wave 
d  =  Drop 

i  =  Initial 

I  =  Interface 

i,0  =  Computation  indices 

L  =  Stagnation  Point 


=  Liquid,  lost,  local 

m        =  Condition  at  maximum 

o        =  Stagnation  condition 

p        =  Point  at  which  pressure  is  acting 

'^^  =  Condition  in  free  stream 


XI 


Abstract  of  Dissertation  Presented  to  the  Graduate  Council 

in  Partial  Fulfillment  of  the  Requirements  for  the 

Degree  of  Doctor  of  Philosophy 


WATER  DROP  DEFORM\TION  AND  FRAGMENTATION 
DUE  TO  SHOCK  WAVE  IMPACT 

By 

William  Edward  Krauss 
June,  1970 

Chairman:   B.  M.  Leadon 

Major  Department:   Aerospace  Engineering 

Results  are  presented  of  experimental  and  theoretical 
studies  of  water  drop  deformation  and  fragmentation  in  the 
high  velocity  air  stream  following  a  normal  shock  wave  moving 
into  still  air.   Experiments  were  conducted  in  a  shock  tube 
with  shock  wave  Mach  numbers  in  the  range  from  1.6  to  3.0  and 
for  drop  diameters  from  .5  to  3.2  millimeters.   From  the  ex- 
perimental data,  it  is  determ.ined  that  the  water  drop  deforms 
axisymmetrically  into  an  approximately  ellipsoidal  shape  which 
is  distorted  as  to  maintain  a  nearly  hemispherical  frontal 
surface  but  with  a  concave  rear  surface.   The  data  also  indi- 
cate that  the  time  for  complete  fragmentation  of  the  water 
drop  is  approximately  one-half  of  the  time  observed  for  com- 
plete reduction  of  the  drop  to  a  trace  of  mist. 

Simultaneous  potential  flow  within  the  distorted  ellip- 
soidal drop  and  mass  loss  by  stripping  of  a  viscous  layer  of 
surface  liquid  are  included  in  a  theoretical  model  for  drop 
deformation  and  fragmentation.   Although  the  actual  drop 
deformation  is  shown  to  deviate  from  the  assumed  ellipsoidal 
shape,  the  general  dimensional  and  displacement  variations 
arc  predicted  by  the  model  and  the  theoretical  mass  history 
of  the  drop  agrees  well  with  the  experimental  data. 
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CHAPTER  I 
INTRODUCTION 

The  deformation  and  fragmentation  of  water  drops  caused 
by  sudden  exposure  to  a  high  velocity  air  stream  currently 
has  important  application  in  the  field  of  missile  and  air- 
craft structural  and  heat  shield  design.   Rain  erosion  by 
the  impingement  of  rain  drops  on  the  exterior  surfaces  of 
missiles  causes  severe  damage  to  ablative  heat  shields  and 
may  result  in  structural  failure  of  the  missile.   The  effects 
of  erosion  damage  can  perhaps  be  minimized  by  judicious 
selection  of  the  heat  shield  material  to  withstand  drop  im- 
pact or  by  suitable  aerodynamic  design  of  the  missile.   In 
this  latter  approach,  the  missile  nose  cap  may  be  designed 
to  generate  a  shock  wave  with  a  standoff  distance  suffi- 
ciently large  to  permit  drop  fragmentation  before  impinge- 
ment upon  the  heat  shield.   Either  of  these  approaches  for 
minimizing  rain  erosion  damage  requires  knowledge  of  rain 
drop  size  and  shape  as  function  of  time. 

Fragmentation  of  liquid  drops  in  high  velocity  gas 
streams  also  occurs  in  fuel  injection  systems.   The  rates 
of  mixing  and  of  combustion  in  such  systems  are  significantly 
enhanced  by  drop  fragmentation. 

A  comprehensive  literature  survey  was  conducted  on 
liquid  drop  deformation  and  fragmentation.   Two  modes  of 
drop  fragmentation  which  occur  are  related  to  Weber's  num- 
ber.  The  first  of  these  modes  occurs  when  the  velocity  of 
the  gas  stream  relative  to  the  drop  is  small  and  corre- 
sponds to  a  magnitude  of  Weber's  number  of  from  1  to  10  or 
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20  times  its  critical  value.   In  this  mode  of  fragmentation, 
the  drop  is  deformed  into  a  "bag"  shape  prior  to  fragmenta- 
tion.  The  second  mode  of  fragmentation  occurs  when  the 
velocity  of  the  gas  stream  relative  to  the  drop  is  high  and 
corresponds  to  the  case  of  Weber's  number  being  very  much 
greater  than  the  critical  value  for  drop  fragmentation.   In 
this  mode  the  drop  is  deformed  axisymmetrically  and  is 
fragmented  by  the  stripping  of  a  liquid  layer  from  its  sur- 
face . 

The  second  mode  of  fragmentation  is  of  interest  for  the 
present  study;  therefore,  only  those  papers  in  the  literature 
which  are  relevant  to  the  second  mode  of  fragmentation  are 
discussed . 

In  general,  the  approach  taken  by  all  investigators  has 
been  similar  in  that  each  used  photographic  techniques  to 
measure  drop  deformation  and  fragmentation  in  the  gas  stream 
behind  a  normal  shock  wave  in  a  shock  tube.   Quantities 
measured  were  drop  diameter,  drop  displacement  and  shock 
wave  Mach  number.   Appropriate  values  of  liquid  density, 
viscosity  and  surface  tension  and  of  gas  density  and  tempera- 
ture were  noted. 

Engel   conducted  comprehensive  studies  wherein  the 
parameters  varied  were  water  drop  size  and  shock  wave  Mach 
number.   These  were  the  first  studies  to  show  that  drop 
fragmentation  is  not  caused  by  the  initial  impact  of  the 
shock  wave  on  the  drop.   In  addition,  these  were  the  first 
studies  to  show  deformation  of  the  drop  into  a  lenticular 
shape  and  stripping  of  a  liquid  surface  layer  from  the  drop. 
Engel  considers  several  phenomena  which  might  explain  the 
stripping  action  but  reaches  no  clear  conclusion  as  to  the 
cause  of  stripping.   In  addition,  an  approach  which  defines 
drop  deformation  based  on  potential  flow  within  the  drop 
is  suggested. 
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Studies  by  Hanson,  et  al.   were  conducted  on  very  small 

drops  (0.5  mm  diameter)  of  various  liquids.  Shock  wave  Mach 
number  was  again  varied  but,  as  in  Engel's  studies,  did  not 
exceed  approximately  1.7.   These  studies  verified  Engel's 

observations  on  largei-  drops. 
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Nicholson  and  Figler   conducted  the  first  studies  with 

shock  wave  Mach  numbers  greater  than  2  thereby  producing  a 
supersonic  air  velocity  behind  the  shock  wave.   This  was 
accomplished  by  using  sub-atmospheric  pressures  in  the 
driven  section  of  the  shock  tube.   The  stripping  phenomenon 
was  again  observed.   Drop  fragmentation  was  shown  to  depend 
largely  on  the  dynamic  pressure  which  was  acting  on  the  drop, 
An  empirical  correlation  for  the  time  required  to  achieve 
complete  drop  fragmentation  is  presented  using  the  data  de- 
veloped for  this  study  as  well  as  those  available  from  Engel 

and  Hanson,  et  al. 

4 
The  studies  of  Ranger  and  Nicholls   were  conducted  on 

water  drops  having  diameters  of  .75  to  4.0  mm.   Shock  wave 

Mach  numbers  from  1.5  to  3.5  in  air  were  used  to  extend  the 

work  of  previous  investigators.   Dynamic  pressures  achieved 

were  significantly  greater  than  those  of  Nicholson  because 

the  pressure  in  the  driven  section  of  the  shock  tube  for 

these  studies  was  one  atmosphere.   The  drop  was  observed  to 

assume  an  ellipsoidal  shape  during  deformation,  and  mass 

removal  was  attributed  to  boundary  layer  stripping.   Using 

empirical  relations  for  drop  diameter  and  drop  relative 

5 
velocity,  the  validity  of  Taylor's   boundary  layer  stripping 

analysis  was  investigated.   Based  on  the  assumption  of  in- 
compressible flow  using  approximate  relations  for  the  ex- 
perimental data,  the  authors  concluded  that  the  agreement 
between  Taylor's  analysis  and  their  experimental  data  was 
encouraging . 


In  the  following  chapters  are  reported  the  results  of 
similar  measurements  which  were  being  started  as  the  Ranger 
and  Nicholls   paper  appeared  in  the  literature.   This  is 
followed  by  an  analytical  treatment  of  simultaneous  deforma- 
tion and  fragmentation  of  the  drop  and  by  a  comparison  of 
the  experimental  results  with  the  analysis. 


CHAPTER  II 

EXPERIMENTAL  INVESTIGATION 

Experimental  Apparatus 

The  experimental  arrangement  used  for  this  investiga- 
tion is  shown  in  Figure  1.   A  stream  of  water  drops  was  in- 
jected by  means  of  hypodermic  needles  into  a  shock  tube  test 
section  which  contained  photographic  windows.   Figures  2  and 
3  show,  respectively,  the  water  drop  injectors  and  the  test 
section.   The  drops  fell  vertically  into  the  test  section 
which  was  open  to  the  atmosphere.   Initial  pressure  in  the 
test  section  was  one  atmosphere.   The  basic  shock  tube  as 
described  by  Clemens   was  modified  by  increasing  the  driver 
length  to  36  inches.   This  was  necessary  to  ensure  that  there 
was  sufficient  time  for  complete  drop  fragmentation  prior  to 
the  arrival  of  the  expansion  wave. 

Both  helium  and  air  were  used  as  driver  gases.   A  model 

600  Honeywell  Strobnar  photographic  flash  unit  was  used  to 

back  light  the  drops  as  is  shown  in  Figure  4.   A  model  ID 

TRW   image  converter  camera  was  employed  to  photograph  the 

drops  after  shock  wave  impact.   The  image  converter  camera 

was  electronically  triggered  by  a  platinum  film  velocity 

g 
sauf^e   located  a  known  distance  from  a  fiducial  marker 
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scribed  on  the  test  section  window.   In  addition,  one  window 
contained  a  millimeter  grid  which  was  used  during  data  re- 
duction to  establish  initial  drop  diameter.   The  experi- 
mental procedure  consisted  of  obtaining  a  history  of  the 
deforiration  and  displacement  of  an  individual  drop  by  taking 
a  series  of  shadow  photographs  at  various  time  intervals 
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after  the  shock  wave  impinged  upon  the  drop.   The  timing 


sequence  was  set  by  a  combination  of  TmV  trigger  delay  gen- 

9 
erators,   Techti-oni: 

shown  in  Appendix  I 
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erators,   Techtronix  oscilloscopes,    and  delay  circuits  as 


Data  Reduction 
The  individual  photographs  were  first  examined  with  a 
Gaertner  optical  micrometer  to  establish  a  history  of  the 
displacement  and  major  diameter  variation  of  the  drop.   The 
time  at  which  the  first  picture  was  taken  was  established  by 
calculating  the  time  required  for  the  shock  wave  to  travel 
from  the  velocity  gauge  to  the  drop  whose  initial  position 
relative  to  the  fiducial  marker  was  recorded  on  the  first 
photograph.   The  time  interval  between  subsequent  photographs 
was  established  by  the  intervals  manually  set  into  the  timing 
circuit.   Following  initial  examination  of  the  photographs 
with  the  optical  micrometer,  the  photographs  were  enlarged 
with  a  Beseler  Vu-Lyte  III  opaque  projector  and  measurements 
were  made  to  determine  the  shape  of  the  windward  surface  of 
the  drop. 

Data  Analysis  and  Results 
The  results  discussed  here  are  for  experiments  that 
cover  the  shock  wave  Mach  number  range  of  M   =  1.6  to  3.0 

o 

in  air  with  drops  having  diameters  in  the  range  D  =  0.5  to 
3.2  millimeters.   Deformation  and  displacement  data  obtained 
from  photographs  were  reduced  in  terms  of  dimensionless  dis- 
placement, X  ,  dimensionless  semi-major  diameter,  b,  dimen- 
sionless  semi-minor  diameter,  a,  and  dimensionless  time,  T. 

Typical  image  converter  camera  photographs  of  shock  wave- 
water  drop  interactions  are  shown  in  Figure  5.   The  dimen- 
sionless and  corresponding  real  times  at  which  the  exposures 
were  made  are  shown  in  Table  I.   Figure  5a,  exposure  5, 
Figure  5b,  exposure  3  and  Figure  5c,  exposure  2  show  the  drops 


at  approximately  their  maximum  deformation  which,  from  Table  I, 
is  seen  to  occur  in  the  time  interval  T  =  1.4  to  1.8.   Figure 
5b  exposure  5  and  Figure  5c  exposure  4  show  similar  stages  of 
drop  fragmentation  and  indicate  that  while  time  T  is  similar, 
actually,  in  real  time,  the  smaller  drop  shatters  much  more 
rapidly  than  the  larger  drop. 

Two  stages  of  development  can  be  observed  in  the  deforma- 
tion and  fragmentation  process.   During  the  first  stage, 
deformation  by  flattening  and  curving  is  the  pi-iraary  phenome- 
non observed  with  little  or  no  fragmentation  in  the  form  of 
mist  visible.   Deformation  into  a  curved  axisymmetrical 
ellipsoid  is  the  result  of  the  external  pressure  distribution 
on  the  drop.   The  second  stage  is  characterized  primarily  by  a 
surface  stripping  action  which  is  due  to  the  shearing  action 
of  the  external  flow  field.   This  shearing  action  rapidly  re- 
duces the  drop  to  a  cloud  of  finely  divided  particles  v^ich 
will  be  termed  "microraist."   The  second  stage  is  well  estab- 
lished by  T  =  .20. 

From  Figure  5  it  is  seen  that  a  well  defined  wake  is 
formed  behind  the  drop  during  early  stages  of  fragmentation. 
The  , shape  of  this  wake  is  similar  to  that  developed  behind  a 
hypersonic  blunt  body  where  the  flow  converges  because  of 
strong  lateral  pressure  gradients  to  form  a  recompression  neck 
region  several  body  diameters  downstream  from  the  rear  stagna- 
tion point.   This  phenomenon  was  also  observed  by  Ranger  and 

4 
Nicholls.    Since  the  liquid  surface  layer  which  is  continu- 
ously being  stripped  from  the  drop  is  able  to  follow  the 
streamline  pattern  of  the  wake,  it  is  reasonable  to  conclude 
that  the  stripped  liquid  is  composed  of  very  fine  water  parti- 
cles, a  micromist.   If  large  pieces  of  liquid  were  being 
stripped  from  the  drop,  the  inertia  of  these  pieces  would  pre- 
clude their  following  the  stream  lines  and  a  significantly 
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different  wake  shape  would  be  observed  in  the  photographs. 
Figure  5c  and  other  similar  photographs  which  shov/  considerable 
portions  of  the  wake  indicate  that  fragmentation  is  a  continu- 
ous process  which  ultimately  reduces  the  entire  drop  to  a 
micromist  cloud. 

Taking  the  time  for  complete  drop  fragmentation  to  be 
that  to  leave  ozily  a  trace  of  mist,  as  was  done  by  Engel   and 
other  investigators,  curves  of  the  deformation  of  the  drop 
and  the  displacement  of  the  drop  were  obtained  and  are  shown 
in  Figures  6  and  7,  respectively,  as  functions  of  dimension- 
less  time,  T.   These  figures  also  include  a  band  which  indi- 
cates the  range  of  experimental  data  available  in  the  literature 
It  is  seen  that  the  data  of  this  study  agree  well  with  those  of 
other  investigators.   Figure  6  indicates  that  complete  frag- 
mentation of  the  drop  occurs  at  approximately  T  =  5.0  as  was 
also  found  by  others. 

Examination  of  the  enlargements  of  the  photographs  ob- 
tained with  the  opaque  projector  indicated  that  significant 
irregularities  occur  in  the  windv/ard  surface  of  the  drop  at 
dimensionless  times  between  approximately  T  =  2.0  and  T  = 
2.5  and  suggests  that  a  breakup  criterion  other  than  reduction 
of  the  drop  to  a  trace  of  mist  might  be  identified.   Measure- 
ments made  on  the  windward  surface  of  the  drop  in  the  various 
photographs  using  the  parameters  as  shown  in  Figure  S  indicate 
this  surface  is  approximately  that  of  an  axisymraetrical 
ellipsoid, the  eccentricity  of  which  varies  with  time  as  shown 
in  Figure  9.   These  data  were  fitted  with  a  mean  fifth  degree 
curve  of  least  mean  square  deviation  expressed  by 

€   =   3.432T-6.515T^  +  4.656T^ 

-1.477T   +  .1753T  (II-l) 

Figure  9  shows  the  eccentricity  varies  from  0  to  a  maxi- 
mum of  .59  and  returns  to  0.   The  frontal  surface  changes  from 
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that  of  a  sphere  to  that  of  an  axisymmotric  ellipsoid  and  re- 
turns to  the  spherical  surface  shape.   The  measured  maximum 
diameter  of  the  drop  and  the  measured  displacement  of  the 
stagnation  point  of  the  drop  were  fitted  with  least  mean  square 
curves  given  by 


and 


B  =  1  +  2.615T-.6776T^-.0971T^  (II-2) 

X   =  .409T^  +  .317  T^  (II-3) 

and  are  shown  in  Figures  10  and  11,  respectively. 

The  fitted  curve  for  the  displacement,  Equation  (II-3), 
also  satisfies  the  known  conditions  at  T  =  0  that  the  dimen- 
sionless  velocity  of  the  drop  is  zero  and  the  dimensionless 
acceleration  is  .818.   The  dimensionless  acceleration  is,  of 
course,  readily  calculated  for  a  sphere  with  a  known  pressure 
distribution.   Details  of  this  calculation  as  well  as  a  method 
for  determining  the  pressure  distribution  on  axisymmetrical 
ellipsoidal  bodies  are  presented  in  Chapter  III. 

The  dimensionless  mass,  m,  of  the  drop  and  the  thickness, 
2a  ,  of  the  drop  were  estimated  by  dividing  the  drop  into 
annular  elements  of  fixed  thickness    y  and  fixed  radius  y 
and  applying  Newton's  Second  Law  to  each  element.   The  dimen- 
sionless force  acting  on  an  element  is  given  by 

F   =   27/ y  Zi.  y  .^p  ,  (II-4) 

where   p  is  the  pressure  acting  on  the  element  at  radius  y  . 
The  dimensionless  mass  of  the  cylindrical  element  is  given  by 

m   =   47ry   -^Y^,  (II-5) 

At  time  T  the  dimensionless  location  of  the  element  with  re- 
spect to  the  fiducial  marker  is  given  by 

X    =   X,   +   -^  (II-6) 

p       L      pp 

which  in  terms  of  the  eccentricity  and  equation  for  an 
ellipse  can  be  written  as 

Xp   =  X^   +  b  (1-£^)"2^1-  [l-(yp/b)^]  2J         (ii_7) 
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Since  X  ,  €  and  b  are  known  functions  of  time,  T,  as  given  by 

i-j 

equations  (II-l),  (II-2),  and  (II-3),  and  since  y   is  known 

for  any  element,  the  accelei-ation  of  the  element  at  any  tine, 

T,  can  be  found  by  differentiating  equation  (II-7)  twice  v/ith 

respect  to  T.   Making  the  substitutions  noted  and  performing 

the  differentiation,  the  diraensionless  acceleration,  X  ,  of 

p 

the  element  is  given  by 

X    =   X,   +   ( bE   +   2bE   +  bE )  ( 1-  5> ) 

+2bE(^^-l)/§  -b^  E(^^-l)/^b         (II-8) 
+bE(^^-l)/S  +   b^E(S>2_i)2/^3^ 


2  1 
where  E   =   (l-€  )2 

§  =£l-(y^/b)2j  2 
Applying  Newton's  Second  Law  and  using  equations  (II-4) ,  (II-5), 
and  (II-7),  the  dimensionless  thickness  of  the  element  is 
given  by 

a^   =   ^p/2   X^  (II-9) 

The  dimensionless  total  mass  of  the  drop  is  readily  found  by 
summing  the  masses  of  the  elements  as  found  by  equation  (II-5). 

Figure  12  presents  the  dimensionless  semi-minor  axis  of 
the  drop,  a  ,  and  the  measured  semi-minor  axis,  a,  of  the  sur- 
face visible  in  the  photographs  which  is  hereafter  called  the 
frontal  surface.   The  water  drop  itself  is  deformed  into  a 
curved  ellipsoid  which  has  the  observed  frontal  surface  but 
a  rearward  surface  which  is  increasingly  concave.   The  axis 
of  the  axial  section  of  this  curved  ellipsoid  is  calculated 

from  the  data.   The  semi-minor  axis  a,  is  calculated  as  indi- 

d 

cated  above,  and  the  semi-major  axis  b   is  measured  along  the 
locus  of  centroids  of  the  cylindrical  segments  whose  lengths 
have  been  calculated  and  plotted  in  Figures  13-24. 

Figure  25  presents  the  semi-major  axis  of  the  drop  as  a 
function  of  dimensionless  time,  T,  and  Figure  26  presents  the 
eccentricity  of  the  drop  as  a  function  of  T  also.   Comparison 
of  Figures  10  and  25  shows  the  semi-major  axis  of  both  the 
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frontal  surface  and  the  drop  reach  maximum  values  at  approxi- 
mately T  =  1.5,   Figure  24  shows  the  cross  section  of  the  for- 
ward portion  of  the  frontal  surface  and  the  deformation  of  the 
water  drop  as  determined  by  the  solution  of  equations  (I 1-4) 
through  (II-9)  for  various  values  of  T.   This  figure  shov/s 
the  water  drop  deforms  into  an  axisymmetric  body  v/hich  is 
crescent  shaped,  or  curved  elliptical,  in  cross  section.   The 
semi-minor  axis  of  the  curved  ellipse  decreases  with  time  as 
mass  is  lost  from  the  drop. 

Figure  27  presents  the  velocity  of  the  stagnation  point 
of  the  drop.   The  velocity  of  the  stagnation  point  is  given 
by 

W   =  .S18T  +  .951T^  (11-10) 

Li 

which  is  directly  obtained  by  differentiation  of  equation 
(II-3). 

Figure  28  presents  the  dimensionless  mass  of  the  drop 
normalized  to  its  initial  mass  and  shows  that  the  mass  of  the 
drop  is  i-educed  to  approximately  2  percent  of  its  initial  mass 
by  time  T  =  2.5  and  indicates  that  for  all  practical  engineer- 
ing purposes  drop  breakup  is  complete  at  approximately  this 
time.   This  result  indicates  the  drop  is  fragmented  in  approxi- 
mately one  half  the  time  for  breakup  as  found  by  using  the 
criterion  that  the  drop  is  reduced  to  only  a  trace  of  mist  as 
observed  on  photographs. 

Experimental  Error 
During  the  experiment,  shock  tube  driver  pressures  were 
measured  within  -3.5  percent  with  a  Heise  model  C-55950  pres- 
sure gauge,  temperatures  within  -1.5  percent  and  velocity  of 
the  shock  wave  within  -3.0  percent.    Time  intervals  between 
photographic  exposures  were  established  by  use  of  the  time 
delay  circuits  in  the  Tektronix  oscilloscopes  which  have  an 
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accuracy  of  -2  percent  of  full  scale  and  the  time  delay  cir- 
cuits in  the  TRW  triga;er  pulse  generators  which  have  an 
accuracy  of  0.01  percent  of  full  scale.^ 

Repeated  measurements  made  with  the  optical  micrometer 
indicated  that  distances  between  two  points,  e.g.,  between 
the  fiducial  marker  and  the  forward  stagnation  point  of  the 
drop,  could  be  measured  within  -3,5  percent.   Comparison  of 
identical  measurements  made  on  the  photographs  using  the 
optical  micrometer  with  those  made  on  the  enlargements  of  the 
photographs  indicate  that  the  latter  measurements  were  accu- 
rate to  within  -7  percent. 

The  standard  deviations  associated  with  each  set  of  ex- 
perimental data  v/hich  were  fitted  with  least  mean  square 
curves  are  shown  in  Figure  9,  10  and  11  and  are  .041,  .085, 
and  .112  for  the  eccentricity,  the  semi-major  diameter  of  the 
front  surface,  b,  and  the  displacement  of  the  stagnation  point, 

X  ,  respectively.   Since  the  standard  deviation  is  the  most 
J-i 

probable  expected  error,  it  is  indicative  of  the  error  due  to 
curve-fitting  the  data.   The  standard  deviation  of  the  dis- 
placement curve,  Figure  11,  and  the  associated  equation, 
(II-3),  are  of  particular  interest  since  equation  (II-3)  is 
differentiated  to  find  the  velocity  and  acceleration  of  the 
stagnation  point  of  the  drop.   The  maximum  standard  deviations 
for  the  velocity  and  the  acceleration  are  calculated  to  be 
.124  and  .146  respectively.   Therefore,  based  on  the  average 
acceleration  of  the  drop  over  the  period  T  =  0  to  T  =  2.5, 
the  calculated  mass  of  the  drop  is  accurate  to  within 
-10  percent. 


CHAPTER  III 

THEORY 

General 
For  convenience  the  analysis  is  divided  into  two  parts. 
In  the  first  part,  the  passage  of  the  shock  wave  over  the 
drop,  the  velocity  and  displacement  of  the  water  drop  are 
mathematically  characterized  assuming  the  flow  remains  at- 
tached to  the  drop  as  the  shock  wave  passes  over  it.   No 
fluid  motion  within  the  drop  is  assumed  to  occur.   This  lat- 
ter assumption  appears  to  be  justified,  based  on  the  photo- 
graphic data  of  Engel.    In  addition,  for  the  conditions  of 
this  study  the  natural  period,  2^  ,  for  steady  oscillation  of 
the  drop  under  surface  tension  as  given  by  Rayleigh, 


V     =   2?r  £2(rn(n-l)(n+2)}  ^ 


is  at  least  three  orders  of  magnitude  greater  than  the  time 


&-' 


for  the  shock  wave  to  pass  over  the  drop. 

The  second  part  of  the  analysis  considers  drop  deforma- 
tion and  fragmentation  and  begins  when  the  shock  wave  reaches 
the  rearmost  point  on  the  surface  of  the  spherical  drop.   The 
flow  field  around  the  drop  and  the  boundary  layer  are  assumed 
to  be  fully  established.   Considering  the  time  required  for 
the  boundary  layer  to  become  established  as  approximately 

that  required  for  boundary  separation  to  take  place,  it  is 

12 
found  from  Schlichting    that 

T      =   .196/3 
sep        r 

The  maximum  dimensionless  separation  time  for  the  conditions 

-2 

of  this  study  is  T      =   1.34  x  10   .   This  is  also  the  time 
^  sep 
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required  for  the  shock  wave  to  move  approximately  a  quarter 
of  a  body  diameter;  thus,  the  assumption  of  a  fully  estab- 
lished boundary  layer  appears  to  be  reasonable. 

Potential  flow  within  the  drop  is  assumed,  and  deforma- 
tion of  the  drop  is  assumed  to  occur  in  a  manner  similar  to 
that  suggested  by  Engel.    The  analysis  here,  however,  is 
assumed  not  to  be  limited  to  small  changes  in  the  drop  radius. 
Stripping  of  a  liquid  layer  from  the  surface  of  the  drop  as 
suggested  by  Taylor   is  also  assumed.   These  two  flov/  mech- 
anisms are  combined  in  a  model  which  undergoes  simultaneous 
drop  deformation  and  mass  loss.   Axial  symmetry  is  assumed 
throughout  the  analysis. 

The  analyses  are  developed  in  terms  of  non-dimensional 
variables  and  reduced  to  numerical  solution  by  the  standard 
forward  marching  technique  used  on  high  speed  electronic 
computers.   The  necessity  for  computer  solution  is  discussed 
as  the  analysis  is  developed. 

Passage  of  Shock  Wave  Over  the  Drop 
The  coordinate  systems  used  to  define  the  position  of 
the  shock  wave,  the  drop  velocity,  and  the  drop  displacement 
are  shown  in  Figure  29.   The  force  acting  on  the  drop  at 
time  T  is: 

f^.Z57rfc^^C^rn^\-K2_)S)^<i^CT)-^2bUt<^<^(j)ld4>LT)  (III-l) 

o 

where  q„ ,  C   ,  K^  and    are  functions  of  T,   K^  is  a 
2    pm   2  2 

parameter  which  governs  the  difference  between  the  pressure 
ahead  of  and  behind  the  shock  wave  as  it  moves  over  the  drop 
and  is  given  by 

K2   =   1.166  (1-Mg2)  p/qgCp^ 
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Since  the  dimensionless  mass  of  the  drop  is 

-7/ 
m   =  —rr- 
6 

the  drop  acceleration  can  be  found  from  Newton's  Second  Lav/ 

provided  the  parameters  under  the  integral  sign  in  equation 

(III-l)  are  known  functions  of  time.   Since  this  is  not  the 

case,  a  forward  marching  numerical  solution  is  employed  wherein 

the  increment  ^4>    is  taken  sufficiently  small  such  that 

q„ ,  C    K^  and  <h    can  be  assumed  with  small  error  to  be  con- 
2 '   pm   2      ^ 

stant  within  the  interval  {Z^  to  <^  +^<^     .      A  substantial  cumu- 
lative error  may  be  incurred  with  this  process,  of  course,  but 
this  is  not  investigated  here  because  the  objective  is  only  to 
estimate  the  drop  deformation  and  displacement  which  occur 
during  the  shock  passage  time.   The  drop  acceleration  is  then 
found  to  be: 


c^~\.SC^^C^yy,    [Cl-K0-3>M>-.ZSSlu'^4>^ol<!p  (III-2) 


J 


4' 


The  drop  velocity  and  displacement  are  found  in  the 
usual  manner  by  integration  with  the  initial  conditions 
OC  =w   =   X    =OatT   =   0.   Using  Zl<p   to  set  the  com- 
pu tat ion  interval,  the  corresponding  time  interval  is: 

^i  Ms 

and  the  velocity  and  displacement  of  the  drop  center  are 
given,  respectively,  by: 

V^=2(°<j  A-T^  -t-l^.J  (III-4) 
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and 


X.=  ^  (W,^  ^^^#^)  ATc  (III-5) 


The  maximum  pressure  coefficient  is  given  by 


^  1.5 


where  for  subsonic  flow  C    is  unity.   C    is  approximated 

pm  pm 

by  a  linear  increase  from  C    =1.0  to  C    =1.8  for 

pm  pm 

.7  <   M   <  1.3  and  by  C    =1.8  for  II      >  1.3.   The  value 
K  pm  R 

of  C    =  1.8  is  typical  of  that  found  experimentally  by 
pm  ^  J       J 

numerous  investigators. 

Drop  Deformation  and  Fragmentation 
The  coordinate  systems  used  in  this  part  of  the  analy- 
sis are  shown  on  Figures  8  and  30,   The  drop  is  assumed  to  be 
an  axisymmetric  ellipsoid  which  has  been  distorted  so  that  the 
front  surface  is  spherical  in  shape.   The  (X,Y)  coordinate 
system  is  fixed  in  space  while  the  z,  ^,  y-coordinate  system 
is  body-fixed  and  ~^  ,  77  are  the  coordinates  fixed  in  the 
boundary  layer.   Assuming  that  potential  flow  exists  inter- 
nally in  the  drop,  let  the  components  of  velocity  of  any 
liquid  particle  in  the  body-fixed  coordinate  system  be 

u   =   -2fJ^  (III-7a) 

V   =   fy  (III-7b) 

w   =   fz  (III-7C) 

where  f  is  an  undetermined  function  of  time  and  is  zero  at 
T  =  0.  Based  on  the  assumption  of  potential  flow,  there 
exists  a  velocity  potential 

<p  -  .sSir^y''^^^-^^^)  (III-8) 
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such  that 

^--2S;5  (III-9a) 

1^  ^5  2  (iii-yc; 


5^^f^  (III-9b) 


Equations  (III-7)  and  (III-9)  clearly  satisfy  the  continuity 

equation  assuming  no  mass ' loss  from  the  drop.   Taking  the  term 

2     2     li 

,5  (u   +  v""  +  w'")  to  be  small   and  integrating  the  equation 

of  motion  the  pressure  is  found  to  be  given  by 

p  =^  +CCt)  (III-IO) 

Applying  equation  (III-IO)  at  the  stagnation  point,  ~^  =   -a  , 
and  at  the  equator,  y  =  b   ,  of  the  drop,  the  following  re- 
lationship is  obtained: 

-f--#^ 5-  (III-ll) 

where  A  P  is  the  pressure  difference  between  the  points 
0  =  77'  and  (^    =      7/^/2    on  the  spherical  surface  of  the  drop. 
All  parameters  in  equation  (III-ll)  are  unknown  functions  of 
time;  hence,  appeal  is  again  made  to  a  numerical  solution 
where  the  computation  interval  is  sufficiently  small  to  assure 
that  the  equation 


is  valid  with  negligible  error.   The  function  f  may  therefoi-e 
be  constructed  as 
J 


ApA~r 


i.=i.C^?^^l. 
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Since  +2fa   is  the  velocity  of  the  point  (-a, 0,0)  on  the  sur- 
d 


face  of  the  axisymmetric  ellipsoid. 


H_^  4.4^  ^  ,  (III-13) 


a  new  value  for  the  semi-minor  axis,  a^,  at  the  end  of  the  com- 
putational interval  is  readily  determined  from: 

The  pressure  difference,  AP  ,  required  for  the  solution 
of  equation  (III-12)  is  given  by: 

Ap-  c^,Cp^0-K3)  (III-15) 

where  q„ ,  C    and  K^  are  functions  of  time.   K   is  defined 
2    pm       o  '> 

and  discussed  later  in  this  chapter. 

From  Figure  30,  the  semi-major  diameter,  b^,  of  the  drop 
is  the  arc  length  along  the  center  line  of  the  drop  and  is 
given  by 


CbA-.5  7rb,..(l-^-|^  -  flf^)  (in-i« 


where 
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The  mass  of  the  drop  at  time,  T,  is  given  by 

IT\;^=m-_^  -vZ^m^  (III-17) 


where  ZX  m .  is  the  mass  loss  from  the  drop  as  determined  by  a 
modified  version  of  Taylor's   boundary  layer  stripping  analy- 
sis which  is  developed  in  Appendix  II.   The  mass  loss,  Am., 
is  given  by 


-  .  0(.ST/  b,_.  S- (i,  ^S,_,  )  AT;^  (III-18) 


The  radius  of  the  frontal  surface  of  the  drop  is  given  by 


b-  ^  ^C-T-Sm/TtO.^:)^^  ^  -       (III-19) 


Equations  (III-14),  (III-16),  (III-17)  and  (III-19)  de- 
fine the  mass  and  approximate  shape  of  the  drop.  Since  the 
shape  of  the  drop  is  assumed  to  be  fixed  during  the  calcula- 
tion interval,  the  acceleration  of  the  center  of  gravity  and 
the  stagnation  point  of  the  drop  is  the  same  and  is  given  by 
equation  (AII-20)  .  Hence,  the  velocity  of  the  stagnation 
point,  W  ,  is  readily  found  by  equation  (III-4)  and  the  dis- 

placement,  X  ,  by  equation  (III-5). 
Li 
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Alternatively,  as  in  the  analysis  of  the  experimental 
data,  the  velocity  and  displacement  of  the  stagnation  point 
can  be  found  by  considering  a  cylindrical  element  of  mass 
which  contains  the  stagnation  point. 

The  force  acting  on  the  cylindrical  element  is  given  by 
(see  Figure  8) 

r  =  Z//^^Ap  •  (III-20) 

The  mass  of  the  element  is 

m  =  4-?r  ^^a_^  (111-21) 

By  Newton's  Second  Law  the  acceleration  of  the  element  is 
<=<=  .5Ap/CL.j  (III-22) 


The  velocity  of  the  stagnation  point,  \'L  ,    is  now  readily 
found  by  equation  (III-4)  and  the  displacement,  X  ,  by  equation 
(III-5) . 

The  pressure  difference,  AP,  in  equation  (III-22)  is  the 
difference  of  the  pressures  acting  on  the  upstream  and  down- 
stream sides  of  the  drop.   The  pressure  on  the  upstream  side 
of  an  axisymmetric  ellipsoidal  drop  is  given  by 


where  it  has  been  assumed  that  the  local  pressure  on  the 
ellipsoid  is  the  same  as  that  found  on  a  sphere  by  applying 
a  tangent  cone  of  constant  vertex  angle  to  both  the  ellipsoid' 


21 

and  the  sphere.   The  local  pressure  on  the  downstream  side  of 
the  drop  is  given  by 

Pa"  Po^■^^^^P'^^3  (III-24) 

K   and  K   are  parameters  which  occur  if  the  relationships  for 
pressure  distribution  on  the  upstream  and  downstream  sides  of 
a  sphere  considering  both  subsonic  and  supersonic  flow  are 
written,  respectively,  in  the  general  forms 

Cp.=  Cpn.Cl-KjSm^^'^   '  (III-24a) 

and 

Cp=  Cp^K^  (III-24b) 

13 

K   for  supersonic  flow  is  unity,    and  for  subsonic  incora- 

1 
pressible  flow   K   =  2.   On  the  rearward  side  of  the  drop  the 

pressure  is  taken  to  be  constant  and  independent  of  <^ 

For  supersonic  flow  K„  =  0,  and  for  subsonic  incompressible 

flow  K   =  -.5.   Transonic  flow  is  assumed  to  occur  between 

M   =  0.7  and  M   =  1.3.   In  this  regime  a  linear  transition 
R  R 

from  the  subsonic  incompressible  to  the  supersonic  values  of 
K^  and  K   is  assumed;  K   and  K   are  thus  given  by 


where  M   is  the  Mach  number  of  the  air  stream  relative  to  the 
drop . 
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The  equations  for  deformation  and  fragmentation  of  the 
drop  are  solved  by  a  forward  marching  numerical  technique  em- 
ploying a  time  interval  sufficiently  small  that  negligible 
error  is  introduced  in  any  given  time  interval  by  assuming 
time  dependent  variables  are  constant  during  that  time  inter- 
val . 

Theoretical  results  and  experimental  data  are  compared 
in  the  next  chapter. 


CHAPTER  IV 
DISCUSSION  OF  RESULTS 

The  numerical  solutions  for  displacement  of  the  drop  in- 
dicate that  for  the  conditions  of  this  study  the  maximum  time 
for  the  shock  wave  to  pass  over  the  drop  is  1.375  /x seconds 
or  a  dimensionless  time  of  .0505,   The  displacement  of  the 
drop  during  this  time  period  is  .009  body  diameters,  which  is 
not  measurable  on  photographs.   The  major  axis  is  calculated 
to  have  changed  less  than  .0003  percent.   Thus  the  assumption 
made  herein  that  the  drop  remains  spherical  during  this  period 
is  justified,  and  the  observations  of  Engel   and  of  Ranger  and 
Nicholls  that  the  drop  exhibits  no  measurable  displacement  or 
distortion  during  shock  passage  are  confirmed. 

For  the  analysis  of  the  deformation  and  fragmentation  of 
the  water  drop  a  relatively  simple  analytical  model  was 
sought  which  would  approximate  the  observed  displacement  of 
the  stagnation  point  of  the  drop,  show  a  variation  of  the  semi- 
major  diameter  of  the  frontal  surface  of  the  drop  similar  to 
that  observed  in  the  experimental  results  and  show  a  change 
in  mass  of  the  drop  approximately  the  same  as  that  found  from 
the  experimental  data.   In  addition,  the  analytical  model 
should  reasonably  determine  the  velocity  of  the  stagnation 
point  of  the  drop  and  the  semi-minor  and  semi-major  axes  or 
the  eccentricity  of  the  drop. 

The  theoretical  model  employed  in  the  analysis  is  based 
on  the  assumptions  of  potential  flow  within  the  drop  and  mass 
flow  from  the  drop  by  the  stripping  of  a  viscous  surface  layer 
of  liquid  at  the  equator  of  the  drop  which  is  taken  to  co- 
incide with  a  point  90  degrees  from  the  stagnation  point 
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on  the  spherical  frontal  surface  of  the  drop. 

Calculated  results  of  the  theory  are  shown  for  the  cases 

of  shock  wave  Mach  numbers,  M  ,  of  1.6  and  3.0  which  cover 

s 

the  range  of  the  experimental  data  of  this  study. 

The  assumption  of  potential  flow  within  the  drop  elimi- 
nates viscous  flow  considerations  from  the  distortion  and 
accounts  for  mass  continuity  during  distortion  of  the  drop 
analysis.   While  this  assumption  appears  to  be  desirable  in 
this  first  attempt  to  model  simultaneously  drop  deformation 
and  fragmentation,  it  is  seen  by  comparison  of  analytical  and 
experimental  results  for  the  semi-minor  diameter,  a  ,  shown  in 
Figures  31  and  12,  respectively,  that  the  model  semi-minor 
diameter  does  not  decrease  initially  as  rapidly  as  observed 
in  experiment.   This  indicates  that  the  assumed  velocity  po- 
tential for  the  liquid  flow  within  the  drop  is  inexact.   As 
will  be  seen,  however,  this  simple  potential  is  quite  adequate 
for  the  principal  purpose  of  the  theory,  that  is,  the  pre- 
diction of  the  mass  history  of  the  drop. 

Since  the  major  diameter,  B,  of  the  frontal  surface  of 
the  drop  is  directly  related  to  the  semi-minor  diameter  as  is 
shown  by  equation  (III-19)  for  mass  of  the  drop,  it  is  ex- 
pected that  the  analytically  determined  major  diameter,  B, 
will  exhibit  a  less  rapid  initial  increase  than  that  observed 
for  B  in  the  experimental  data.   Comparison  of  Figures  10  and 
32  show  this  to  be  the  case.   This  comparison  also  indicates 
that  the  maximum  analytically  determined  values  of  the  major 
diameter,  B,  are  within  10  percent  of  the  maximum  experi- 
mentally measured  values.   Further,  the  dimensionless  time 
at  which  the  calculated  major  diameter,  B,  reaches  a  maximum 
is  within  15  percent  of  the  experimentally  observed  time  at 
which  B  reaches  a  maximum.   The  differences  between  experi- 
mentally observed  and  theoretical  values  of  a^  and  B  during 
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latter  phases  of  fragmentation  of  the  drop  are  attributed  to 
assumptions  made  in  the  boundary  layer  stripping  analysis  and 
are  discussed  later  in  this  chapter. 

Figure  33  presents  the  calculated  mass  of  the  drop  nor- 
malized to  the  initial  mass  of  the  drop.   Comparison  of  this 
figure  with  Figure  28  for  the  experimentally  determined  nor- 
malized mass  of  the  drop  indicates  good  agreement  between  the 
observed  and  calculated  values.   From  these  figures  it  is  also 
seen  that  the  calculated  mass  of  the  drop  is  somev/hat  less 
than  the  experimentally  observed  value  when  the  mass  of  the 
drop  is  reduced  to  approximately  20  percent  of  its  original 
mass.   From  Figure  31  it  is  seen  that  at  approximately  T  = 
1,7  the  drop  is  becoming  very  thin  and  clearly  the  viscous 
boundary  layer  is  no  longer  small  compared  to  the  thickness 
of  the  drop.   Hence  the  theoretical  model  which  is  based  on 
the  assumption  of  a  thin  boundary  layer  is  inaccurate.   This 
deterioration  of  the  model  for  mass  loss  is  also  reflected 
through  equation  (III-19)  in  the  behavior  of  the  major  diam- 
eter of  the  frontal  surface  as  shown  in  Figure  32  for 
T  >•  1.8  to  2.0.   Similarly,  the  semi-minor  diameter,  a  ,  is 
also  affected  through  equations  (III-16)  and  (III-12) . 

Figure  34  presents  the  theoretical  displacement  of  the 
stagnation  point  of  the  drop.   Comparison  of  Figure  34  and  11 
indicates  that  the  theoretical  model  predicts  smaller  dis- 
placements during  the  breakup  process  than  aie  observed  ex- 
perimentally,   At  T  =  2.5,  however,  the  experimentally 
observed  displacement  falls  between  the  calculated  values 
for  the  range  of  shock  wave  Mach  number  of  this  study.   The 
difference  between  the  experimentally  observed  and  the  cal- 
culated displacement  is  attributed  primarily  to  the  assump- 
tion made  in  the  analysis  that  the  maximum  pressure  coeffi- 
cient varies  linearly  with  the  relative  Mach  number,  M  ,  in 
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the  transonic  regime  and  is  constant  in  the  subsonic  and  super- 
sonic regimes  as  given  by  equation  (III-6). 

This  or  a  similar  assumption  as  to  the  variation  of  the 
local  pressure  coefficient,  C  ,  on  a  sphere  was  necessary 
since  no  experimental  data  on  pressure  distributions  around 
spheres  or  hemispheres  in  the  transonic  regime  appear   to  be 
available. 

Figure  35  presents  the  theoretical  velocity  of  the  stag- 
nation point  of  the  drop.   Since  the  velocity  is  the  time  rate 
of  change  of  the  displacement,  it  is  to  be  expected  from  the 
slope  of  the  displacement  curves  in  Figure  34  that  the  pre- 
dicted velocity  will  be  larger  than  that  experimentally  ob- 
served.  Comparison  of  Figures  27  and  35  show  this  to  be  the 
case  with  the  calculated  velocity  being  approximately  25  per- 
cent greater  than  the  experimentally  determined  velocity  at 
T  =  2.5.   This  difference  is  also  attributed  to  the  assumed 
linear  variation  of  the  maximum  pressure  coefficient  as  dis- 
cussed in  the  preceding  paragraph. 

Figure  36  presents  the  calculated  time  variation  of  the 
eccentricity  of  the  drop.   Comparison  of  this  figure  with 
Figure  14  indicates  similar  timewise  variation  of  the  eccen- 
tricity.  While  both  the  theoretically  predicted  and  experi- 
mentally determined  eccentricities  approach  unity  between 
T  =  .6  and  T  =  .7,  the  experimentally  determined  value  in- 
creases more  rapidly  than  does  the  analytically  predicted 
value.   This  is  due  to  the  relatively  slower  decrease  in  the 
analytically  predicted  variation  in  the  semi-minor  diameter, 
a  ,  of  the  drop  than  is  observed  in  the  experimental  results 
as  was  discussed  earlier  in  this  chapter. 


CHAPTER  V 
CONCLUSIONS 

This  investigation  of  deformation  and  fragmentation  of 
a  water  drop  due  to  shock  wave  impact  has  produced  several 
new  findings  from  both  the  experimental  and  theoretical  re- 
sults . 

The  primary  experimental  data  which  are  measurements  of 
the  displacement  and  the  major  diameter  of  the  deformed  water 
drop  are  in  good  agreement  with  those  of  other  investigators. 

Analysis  of  the  experimental  data  shows  that  the  water 
drop  is  deformed  into  an  axisymmetric  shape  which  has  a 
nearly  spherical  ellipsoidal  frontal,  or  windward,  surface 
and  a  concaved  rear  surface. 

Determination  of  the  mass  of  the  drops  as  a  function  of 
time  from  the  experimental  data  permits  a  new  criterion  to 
be  used  to  establish  the  time  for  complete  fragmentation  or 
breakup  of  the  drop.   Using  the  time  at  which  the  mass  of 
the  drop  is  reduced  to  approximately  1  percent  of  its  origi- 
nal mass,  it  is  found  that  drop  breakup  occurs  at  approxi- 
mately T  =  2.5.   This  is  approximately  one-half  of  the  time 
identified  by  other  investigators  who  used  the  criterion  that 
breakup  time  is  that  at  which  the  water  drop  is  observed  to 
be  reduced  to  a  trace  of  mist  on  photographs  of  water  drop 
fragmentation.   Using  this  criterion  the  time  for  complete 
breakup  of  the  water  is  typically  given  as  T  =  5.0. 

The  approximate  analytical  model  developed  in  this 
study  indicates  simultaneous  deformation  and  mass  loss  from 
a  water  drop  can  be  treated  using  a  simple  theoretical  model. 
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It  is  concluded  from  a  comparison  of  theory  and  experiment 
that  the  classical  Taylor's  boundary  layer  stripping  analysis 
must  be  modified  in  order  to  obtain  a  realistic  prediction  of 
the  mass  history  of  the  water  drop.   The  modification  proposed 
herein  (Appendix  II)  permits  the  viscous  profile  to  extend 
through  the  entire  drop  as  it  becomes  thin. 

It  is  further  concluded  that  while  the  assumption  of 
potential  flow  within  the  drop  is  satisfactory  for  an  approxi- 
mate theoretical  model  as  was  used  herein,  a  more  exact  theory 
of  the  internal  flow  should  be  developed. 

It  is  also  concluded  that  the  general  lack  of  experi- 
mental data  on  the  pressure  distributions  on  spheres,  hemi- 
spheres and  ellipsoids  in  the  transonic  and  low  supersonic 
flow  regime  will  restrict  the  accuracy  of  calculations  on  the 
displacement,  deformation  and  fragmentation  of  water  drop 
which  are  imparted  with  shock  waves  moving  in  the  Mach  number 
range  of  approximately  1.6  to  4.0. 

Finally  it  is  concluded  from  the  good  agreement  between 
the  experimentally  and  theoretically  determined  mass  of  the 
drop  that  the  mass  calculation  is  relatively  insensitive  to 
the  precise  form  of  the  potential  flow  function  used  to  cal- 
culate deformation  of  the  drop. 
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TABLE  I 
TIME  SEQUENCE  FOR  PHOTOGRAPHS  SHOWN  IN  FIGURE  5 


Figure  5a  -  Exposure;  D  =  2.06  mm 
1 
2 
3 
4 
5 


Dimensionless 

Real 

Time,  t, 

Time,  T 

-7 

S( 

3conds 

1 6  mm . 

0 

0 

.460 

44 

.985 

94 

1.140 

109 

1.455 

139 

Figure 

5b  - 

-  Exposure; 

D  = 

=  2.74  mm. 

1 

0 

2 

.220 

3 

1.155 

4 

2.060 

5 

2.980 

Figure 

5c  - 

-  Exposure; 

D  = 

=  .985  mm. 

1 

0 

2 

1.200 

3 

1.850 

4 

2.740 

5 

4.260 

0 

30 
124 
224 
324 

0 

55 

85 

125 

195 
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FIGURE  6,   EXPERIMENTAL  DROP  DEFORMATION  DATA 
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FIGUKE    7.       EXPERIMENTAL   DROP    STAGNATION   POINT 
DISPLACEMENT    DATA 
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FIGURE    8.       COORDINATE    SYSTEMS    FOR  ANALYSIS 
OF   EXPERIMENTAL   DATA 
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FIGURE    10.       MEASURED   DEFORMATION   OF   THE   DROP   AND 

LEAST   MEAN   SQUARE   CURVE   FIT   TO   THE   DATA 
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FIGURE    11 
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MEASURED   DISPLACEMENT    OF   THE    STAGNATION    POINT 
OF   THE   DROP   AND   LEAST    MEAN   SQUARE    CURVE    FIT 
TO   THE   DATA 
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FIGURE    12 
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MEASURED  SEMI-MINOR  DIAMETERS  OF  THE  FRONTAL 
SURFACE  OF  THE  DROP  AND  OF  THE  DROP 
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SCALE:  I- 


H  .lb 


FIGURE  13. 


MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE 
FRONTAL  SURFACE  AND  CALCULATED  CROSS 
SECTIONAL  SHAPE  OF  THE  DROP  AT  T  =  .4 
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FIGURE  14 


MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE 
FRONTAL  SURFACE  AND  CALCULATED  CROSS 
SECTIONAL  SHAPE  OF  THE  DROP  AT  T  =  .6 
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SCALE:  f 1  .lb 


'     1 

FIGURE  15.   MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE 

FRONTAL  SURFACE  AND  CALCULATED  CROSS  SECTIONAL 
SHAPE  OF  THE  DROP  AT  T  =  .8 


45 


SCALE:  I 1  .lb 


FIGURE  16.   MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE 
FRONTAL  SURFACE  AND  CALCULATED  CROSS 
SECTIONAL  SHAPE  OF  THE  DROP  AT  T  =  1.0 
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SCALE:  I 1  .lb 


FIGURE  17. 


MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE  FRONTAL 
SURFACE  AND  CALCLXATED  CROSS  SECTIONAL  SIL-^PE 
OF  THE  DROP  AT  T  =  1.2 
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SCALE: 
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FIGURE  18.   MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE  FRONTAL 
SURFACE  AND  CALCULATED  CROSS  SECTIONAL  SHAPE 
OF  THE  DROP  AT  T  =  1.4 
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SCALE:  t- 
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FIGURE  19.   MEASURED  CROSS  SECTIONAL  SIlAPE  OF  THE  FRONTAL 
SURFACE  AND  CALCULATED  CROSS  SECTIONAL  SHAPE 
OF  THE  DROP  AT  T  =  1.6 
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SCALE:  h 
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FIGURE  20 


MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE  FRONTAL 
SURFACE  AND  CALCULATED  CROSS  SECTIONAL  SHAPE 
OF  THE  DROP  AT  T  =  1.8 
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SCALE:  I 1. lb 


FIGURE  21 


MEASURED  CROSS  SECTIONAL  SHAPE  OF  THE  FRONTAL 
SURFACE  AND  CALCULATED  CROSS  SECTIONAL  SHAPE 
OF  THE  DROP  AT  T  =  2.0 
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FIGURE  22. 


MEASURED  CROSS  SECTIONAL  SIL-iPE  OF  THE 
FRONTAL  SURFACE  AND  CALCULATED  CROSS 
SECTIONAL  SHAPE  OF  THE  DROP  AT  T  =  2.2 
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SCALE: 
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FIGURE   23. 


MEASURED  CROSS  SECTIONAL  SHAPE  OF 
THE  FRONTAL  SURFACE  AND  CALCULATED 
CROSS  SECTIONAL  SHAPE  OF  THE  DROP 
AT  T  =  2.4 
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FIGURE  25.   MEASURED  SEMI-MAJOR  DIAMETER  OF  THE  DROP 
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FIGURE    26. 
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FIGURE    27.       EXPERIMENTALLY   DETERMINED   VELOCITY    OF   THE 
STAGNATION   POINT   OF    THE   DROP 
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FIGURE  28.   EXPERIMENTALLY  DETERMINED  NORMALIZED  MASS 
OF  THE  DROP 
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FIGURE  31.   THEORETICAL  SEMI-MINOR  DIAMETER  OF  TIIE  DROP 
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FIGURE    32.       THEORETICAL   MAJOR  DIAMETER   OF   THE   FRONTAL 
SURFACE   OF   THE   DROP 
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FIGURE  33.   THEORETICAL  NORMALIZED  MASS  OF  THE  DROP 
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FIGURE  34.   THEORETICAL  DISPLACEMENT  OF  THE  STAGNATION 
POINT  OF  THE  DROP 
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FIGURE  35.   THEORETICAL  VELOCITY  OF  THE  STAGNATION 
POINT  OF  THE  DROP 


65 


C3 
C3 


Mj  =  3.0 
Ms  =  1.6 


1.0        1.5         2.0 
DIMENSIONLESS  TIME,  T  -  t,3  U2/D 


2.5 


FIGURE  36.   THEORETICAL  ECCENTRICITY  OF  THE  DROP 
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APPENDIX    II 


BOUNDARY  LAYER  STRIPPING  ANALYSIS 

The  analysis  for  mass  loss  from  the  boundary  layer  on  a 

5 
spherical  drop  is  attributable  to  Taylor.    The  analysis  pre- 
sented here  follows  Taylor  but  is  more  general  in  that  it 
applies  to  axisymmetric  bodies. 

Upon  sudden  exposure  of  a  water  drop  to  a  high  speed  air 
stream  boundary  layers  will  form  both  in  the  air  and  in  the 
water.   The  coordinate  system  used  is  shown  in  Figure  30  where 
■^  is  the  curvilinear  coordinate  along  the  interface  separat- 
ing the  two  fluids  and  >;?  is  the  coordinate  perpendicular  to 
it  and  is  taken  positive  in  either  direction  from  the  inter- 
face.  An  approximate  solution  to  this  double  boundary  layer 
problem  is  obtained  by  assuming  arbitrary  simple  velocity  dis- 
tributions containing  several  unknown  parameters  which  are 
determined  through  use  of  Karman's  momentum  integral  relations 
Assume  that  the  flow  is  steady  and  incompressible  and  that 
for  the  gas 

^--   ^-U'\J.)^ypl-^f^-}  (Aii-i) 


-x 


and  for  the  water 


where,  U.     is  the  local  velocity  in  the  gas  outside  of  the 
boundary  layer,  U   is  a  function  of  C  and  <=>(   and  c<j,  are 
constants.   The  steady  state  boundary  layer  momentum  integral 
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equations   are    for    the   gas 
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NT  1  e^ 


cX  =  o 


r- 


(AII-3) 


and   for   the  water 


M'<''?'*v|r«'^r : 


''O 


d 


r 


(AII-4) 


where  the  pressure  gradient  in  the  liquid  layer  is  given  by 


^? 


? 


(AII-5) 


At  the  gas-water  interface  the  shear  stresses  must  be  equal 
and  therefore 


rA§\.^^MtK-^ 


^y  ^l- 


c)>7  \ 


(AII-6) 
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Applying  equations  (AII-1),  (AII-2) ,  and.  (AII-6)  and  evaluat- 
ing the  result  at  9^  =  ^/2  where  for  a  sphere 

equations  (AII-3)  and  (AII-4)  become,  respectively, 


(j^ylH.^:    _Xy_,a-U.)A       __o  (AII-7, 


Afj'  =<^  ry 


and 


4^  ^^  Vj^ 


where  .^  is  the  kinematic  viscosity  coefficient. 
These  equations  yield  directly 


% 


^2. 


-{4:^/Ur^jz]-  (Aii-S) 


and 


/  '/7 

"<^"=[^  V(^x+-')^^}  (AII-10) 
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Differentiating  equations  (AII-1)  and  (AII-2)  with  re- 
spect to  >?  ,  inserting  the  result  into  equation  (AII-6)  and 
evaluating  at  -^^  =0,  the  following  equation  is  obtained: 

i--x;^--uj4-xj,^(i-p^^^:')-o  (Aii-ii) 


Since  U   is  the  interface  velocity  and  is  expected  to  be  very 
small,  equation  (AII-11)  can  be  approximated  by: 


l-^V^Kr  ''O  '"'''" 


and  the  interface  velocity  is 


^i-f^^'«^r^  ^^''-''^ 


The  mass  of  water  in  the  surface  layer  being  swept  along 
by  the  gas  stream  at  a  distance  equal  to  the  arc  length  be- 
tween the  stagnation  point  and  the  equator  is  given  by 


where  u^  is  taken  to  be  the  local  velocity  as  given  by  equation 
(AII-2)  modified  by  a  drop  distortion  velocity,  U^,  due  to 
internal  flow  within  the  drop.   Hence 
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-r 


u^-^,^^x"^.  ^^p[,;r^j 


(AII-15) 


The  distortion  velocity  U   is  assumed  to  be  proportional 

to  the  velocity  of  the  fluid  in  the  drop  at  the  equatorial 

point  such  that 

U,   =   .065v 
d 

where  v  is  given  by  equation  (III-7b)  and  .065  is  an  em- 
pirical correlation  factor.   Equation  (AII-14)  can  be  written 
as 


CO 


-^/LU't 


.0(o5 


-h  fu^U.  exp>— ^-TCT? 


^o 


'X^o 


(AII~16) 


where  the  first  integral  is  taken  over  the  boundary  layer 
thickness,  S  ,  defined  such  that 


exp 


r-S 


i<=<o  f^  3 


=  .01 


^^? 


Performing  the  integrations  indicated  in  equation 
(AII-16)  there  results 


^  -abloksiv-i  H^)\^.--^} 


(AII-17) 
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Integrating  equation  (AII-17)  over  a  small  time  interval, 
A  T,  so  that  all  time  dependent  parameters  can  be  considered 
constant  with  negligible  error  and  applying  equations  (AII-9) 
and  (AII-13)  to  the  result  the  mass  loss  from  the  drop  is 
found  to  be 


-  •^^5?^^d-.S.(i+i:.J^T^        '    '(Aii-18) 


The  local  velocity,  IL,  ,  at  the  equator  of  the  drop  is 
approximated  by 


"^  =  K^  (I-  P^J//3)  (AII-19) 


Equation  (AII-19)  is  based  on  the  assumption  of  isentropic 

flow  around  the  drop.   This  assumption  is  of  course  not 

strictly  valid  for  supersonic  flow;  however,  since  M  ■^  1.35 

R 

for  the  present  study,  the  entropy  change  is  small  and  the 
assumption  of  isentropic  flow  introduces  little  error.   For' 

subsonic  incompressible  flow,  M   <  0.7,  the  velocity  at  the 

14  ^ 

equator  of  the  drop    is 


U,=L5(l-(DV/y(S 


and  for  supersonic  flow  with  M   >  1.3 

R 


t^-Cl-<^v</)/(3 
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In  the  Mach  number  range,  .7  <  M  -^  1.3  K   is  assumed 
to  vary  linearly  and  is  given  by 

K  =  1.1(2  5 -O.ISrsCi'Mi^-.?!-  ll^R-1-3 


The  velocity,  W.,  appearing  in  equations  (AII-18)  and 
(AII-19)  is  the  velocity  of  the  drop  assuming  no  deformation 
is  occurring  and  is  given  by  equation  (III-4)  with: 
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